Intervalley scattering and localization behaviors of spin-valley coupled Dirac fermions 
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We study the quantum diffusive transport of multi-valley massive Dirac cones coupled by in- 
tervalley spin-orbit scattering. We show that the intervalley spin-orbit scattering and intravalley 
spin-conserved scattering can be distinguished from the quantum conductivity that corrects the 
semiclassical Drude conductivity, due to their distinct symmetries and localization trends. In im- 
mediate practice, it allows transport measurements to estimate the intervalley scattering rate in 
hole-doped monolayers of group- VI transition metal dichalcogenides (e.g., molybdenum dichalco- 
genides and tungsten dichalcogenides), an ideal class of materials for valleytronics applications. The 
results can be generalized to a large class of multi- valley massive Dirac systems where time-reversal 
symmetry demands opposite spins in opposite valleys. 

PACS numbers: 73.63.-b, 75.70.Tj, 85.75.-d 



Valley index associated the degenerate energy extrema 
of Bloch bands in momentum space is a well defined de- 
gree of freedom for low-energy carriers which can encode 
information just like the spin index in spintronics. This 
has led to the concept of valleytronics, a new type of elec- 
tronics based on manipulating the valley index of carri- 
ers [THS]. The timescale of intervalley scattering deter- 
mines how long the information represented by valley po- 
larization can be retained. For valleytronics applications 
to be practical, this timescale shall be long as compared 
to the typical time for the control of valley dynamics. 
Measurement of intervalley scattering time is thus criti- 
cal for identifying potential valleytronics materials. 

Several extensively studied monolayer 2D crystals are 
promising materials to host valley-based electronics, in- 
cluding graphene, and graphene-like crystals such as sil- 
icene [6T-HU] and monolayer group- VI transition metal 
dichalcogenides [TTHH] . In these hexagonal 2D crystals, 
both the conduction and valance band edges are at the 
two incquivalcnt valleys at K points (corners of first Bril- 
louin zone), which are related by time reversal. In mono- 
layer dichalcogenides, ground breaking theoretical and 
experimental progresses on dynamical control of valley 
are recently achieved [T5HTS] . Unlike graphene, mono- 
layer dichalcogenides are described by massive Dirac 
fermions, and intrinsic spin-orbit coupling (SOC) gives 
rise to a spin splitting which must have opposite sign 
at the two valleys by the requirement of time reversal 
symmetry (see e.g. Fig. 1). This effective coupling be- 
tween the spin and valley indices can have two significant 
consequences for valleytronics: (i) the interplay between 
spin and valley degrees of freedom; (ii) a unique form of 
intervalley scattering which must be accompanied by a 
simultaneous spin flip. 

In this Letter, we study the quantum diffusive trans- 
port of multi-valley massive Dirac fermions with spin- 
valley coupling. Without loss of generality we choose the 
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FIG. 1: Low-energy effective band structures of MX2 mono- 
layer in K and K' valleys. The Fermi surface (at the top 
edges of shadowed areas) intersects the highest two valance 
bands of opposite spin orientations (marked by f and J.). Hor- 
izontal dotted arrows guide the intervalley spin-flip scattering 
assisted by spin-orbit interaction. 



model of monolayer dichalcogenides MX2 (M = Mo, W; 
X = S, Se) for a concrete discussion. This problem is the- 
oretically unique from the conventional 2D electron gas 
systems or the massless graphene systems without spin- 
orbit coupling. We show that the spin-valley coupled 
band edge results in distinct symmetries and localization 
behaviors in intra- and inter-valley scattering dominant 
regimes. Namely, the spin-conserving intravalley scatter- 
ing leads to negative quantum conductivity and positive 
magnetoconductivity from weak localization (WL) effect, 
whereas the intervalley scattering assisted by spin-orbit 
interaction gives rise to positive quantum conductivity 
and negative magnetoconductivity from weak antilocal- 
ization (WAL). The quantum conductivity and magne- 
toconductivity have logarithmic dependence on the ratio 
between intra- and inter-valley scattering rates. Such 
a dependence makes possible measurement of interval- 
ley scattering rate in a vastly broad range in monolayer 
dichalcogenides, from three orders of magnitude below to 
three orders of magnitude above the intravalley scatter- 
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ing rate, which can be readily tested in available dichalco- 
genide field effect transistors [IMSTj . Experimental veri- 
fications of these phenomena will reveal the nature of im- 
purities and provide general guidance on the suppression 
of intervalley scattering for a longer valley life time nec- 
essary for valleytronic applications. Our approach and 
results can also be readily extended to spin-valley cou- 
pled Dirac fermions with much smaller masses in silicone 
or monolayers of germanium and tin [5]. 

The band edge electrons and holes in monolayer 
MX2 are well described by the massive Dirac fermion 
model [15] : 

H = hv{±k x a x + kydy) + ~a z ± Xs z ® ° z , (1) 

where ± stands for K and K' valleys, respectively. The 
Pauli matrices cr x ,y,z act on the pseudospin indexing the 
A = d z 2 and B = (d x 2_ y 2 ± id xy )/^/2 orbitals while 
s z =t,4- is the z-component of real spins. (k x ,k y ) is the 
wave vector measured from K (K 1 ) points. The Hamilto- 
nian is described by three parameters: the SOC strength 
A, the band gap A, and the effective velocity v. The 
resulting band structure consists of four sets of massive 
Dirac cones, two in K and two in K' valleys. An impor- 
tant feature is the large SOC splitting [15l|22](2A - 0.15 
cV in molybdenum dichalcogenides and 2A ~ 0.4 eV in 
tungsten dichalcogenides) between the spin-up and spin- 
down states at the valence band top. The conduction 
band bottoms remain degenerate. The band dispersion 
is schematically shown in Fig. [I] 

Because of the large spin splitting in the valence bands, 
in the following we focus on the localization effect in hole- 
doped samples. It is instructive to first consider two lim- 
iting regimes when either intra- or inter-valley scattering 
dominates. 

(1) Intravalley scattering dominant regime. In the va- 
lence bands, when holes are confined to a single valley, 
the SOC acts like a strong Zeeman field, which fully po- 
larizes the spins. Therefore, the intravalley scattering 
can only happen between the same spin species. At the 
same time, because of the large band gap (A ~ 1.6 eV), 
the orbital pseudospin is almost fully-polarized and does 
not play a role. The frozen pseudospin can be quanti- 
tatively seen from the Berry phase [23], which for the 
valence bands is given by 
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where V measures the Fermi energy from the valence 
band top. A is much larger than A and allowed V 
(e [0, 2A] eV) in MX2 [15] . The spin-conserved intraval- 
ley scattering that also conserves time-reversal symmetry 
belong to the orthogonal symmetry class [24] , which gives 
rise to WL [25]. 

(2) Intervalley scattering dominant regime. The in- 
tervalley scattering in MX2 must break spin-rotational 
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FIG. 2: (a) Sketch of band structure, where V £ [0, 2A] eV 
measures the Fermi energy from the valence band top. [(b)- 
(d)] Magnetoconductivity Act under a perpendicular magnetic 
field B z , and when the Fermi energy intersects the two higher 
valence bands. The positive (negative) logarithmic magne- 
toconductivity is the signature of the weak localization (an- 
tilocalization). (b) No intervalley scattering (tq/ti = 0) for 
different V. (c) At V = 0.1 eV and for different to/tt, the 
ratio between intervalley and intravalley scattering rates, (d) 
Zoom-in of (c) for tq/tj = 0.1,1,10. Parameters: A — 1.66 
eV and A = 0.075 eV [15], mean free path £ = 10 nm, and 
phase coherence length i^, = 300 nm. 



symmetry because of the opposite spin splitting at the 
two valleys (Fig. [I]). For two possible spin-flip mecha- 
nisms, the magnetic scattering and spin-orbit scattering, 
the latter is likely to exist in the materials. The spin-orbit 
scattering preserves time-reversal symmetry, along with 
the broken spin-rotational symmetry, it has the symplec- 
tic symmetry [23] which leads to WAL [25] , 

To test the above picture, we calculated the WL and 
WAL using the diagrammatic techniques 26 30j . Be- 
cause the results for gapless graphene without spin-orbit 
coupling can not be applied directly, we generalize the 
calculation to this multi-valley massive Dirac cones cou- 
pled by intervalley spin-orbit scattering. The method is 
based on the perturbation expansion of the Kubo formula 
for conductivity. The Hamiltonian in Eq. ([I]) is taken as 
the unperturbed part and scattering potentials as pertur- 
bations. The scattering potentials are given as follows. 
The spin-conserved intravalley scattering is modeled by 
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where u is the potential of an impurity at position R. 
The subscript k, k' means the potential scatters an elec- 
tron of wavevector k' to k. The spin-flip intervalley scat- 
tering is given by spin-orbit scattering 
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U A ' B = i J2 «V (k '- k) ' R k xk'-z, 

ReA/B 

U£ /B = i « R e 4(k '- k) ' R k x k' ■ (x±iy) , (4) 

RGA/B 

where the matrix is spanned by the basis functions in the 
order of |+t^4), |+t-B), \~l A )i \~\-B), and x,y, z stand 
for the unit vectors. The extra i in U^ s £i protects time- 
reversal symmetry in spin-orbit scattering as i, k, k', and 
spin change sign under time reversal. For mathematical 
simplicity, we assumed short-range and delta correlation 
for both U° and U l ' 3 °. Although intravalley scattering 
should be related to long-range potential, the practice 
by the delta potential was justified [31]. 

From the diagrammatic techniques, the total conduc- 
tivity is given by 

a = * D + * F , (5) 

where a D is the semiclassical (Drude) conductivity. Ac- 
cording to the Einstein relation, a D is given by 

a D = e 2 NpD (6) 

where Np — Ep/(-Kh 2 v 2 ) is the density of states at the 
Fermi surface with Ep the Fermi energy measured from 
the Dirac point, D = t] v v f t/2 is the diffusion constant 
with rj v the ladder diagram correction to the velocity of 
Dirac fermions [32 , vp the Fermi velocity, and the total 
scattering time given by r = (1/t + 1/ti) , where r 
and ti are intravalley and intervalley scattering times, re- 
spectively. Therefore, a D is a function of the total scat- 
tering time, thus can not distinguish the contributions 
from intervalley and intravalley scattering. Besides, a D 
is insensitive to magnetic field. For these reasons, the 
quantum interference correction to the Drude conductiv- 
ity a F (quantum conductivity for short) has to be con- 
sidered, which is found to take the form of logrithmic 
functions 

nh y AV 2 + £~ 2 Xf 2 + I- 2 J W 

where C\ and Co are weight factors for the intravalley and 
intervalley contributions. Xq and X\ are corresponding 
characteristic lengthes that effectively reduce the mean 
free path £ and phase coherence length l^. The inter- 
valley/intravalley scattering ratio is incorporated in the 
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FIG. 3: Zero-field quantum conductivity a F as a function 
of temperature T for different to/ti. Positive (negative) a F 
corresponds to weak antilocalization (localization). Other pa- 
rameters are the same as those in Figs. [2](c) and (d). 

analytic expressions for Cu and Xo,i [33] - Due to its 
quantum interference origin, a F can be suppressed by a 
perpendicular magnetic field B z , leading to a magneto- 
conductivity 

Aa(B z ) = a F (B z ) - a F (0) = ^r(C F + M), (8) 

irn 

where = ^(£ 2 B /£ 2 + ^yH£ 2 B /£f), l/ lj = 1/%+ 1/Xf, 

\E' is the digamma function, and £b = \/^-/(4e|£^|) is the 
magnetic length [55] . 

Figure [2jb) shows the magnetoconductivity in the ab- 
sence of intervalley scattering, and when the Fermi en- 
ergy intersects with the two highest valence bands. The 
WL is a quantum interference induced suppression of con- 
ductivity, which can be lifted by magnetic field and gives 
positive magnetoconductivity. Figure [2jb) shows that 
the magnetoconductivity is always positive in the allowed 
range of the Fermi energy, indicating the WL from each 
valence band. Figures [2]Jc) and (d) shows the magneto- 
conductivity as the intervalley scattering increases. It is 
convenient to define the ratio between intervalley spin- 
orbit and intravalley scattering rates 

ro/n, (9) 

which increases with increasing intervalley scattering. In 
Figs. |2jc) and (d), as tq/t\ increases, the magnetocon- 
ductivity changes from positive to negative, correspond- 
ing to a crossover from WL to WAL. The magnetocon- 
ductivity covers a wide range of tq/ti, from < 0.001 to 
> 1000. When t /ti > 10 or < 0.1, the magnetocon- 
ductivity changes sensitively with Tq/ti. This sensitivity 
allows us to estimate even a very small t\ when To dom- 
inates, or vice versa. Considering the low mobility in 
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FIG. 4: Weak localization (WL) and antilocalization (WAL) 
in the conduction bands. Spin-up (t) and spin-down (J.) are 
offset for clarity. Solid and dotted arrows represent spin-orbit 
and spin-conserved scattering, respectively. 



MX2 [20] , the calculation has assumed a short total scat- 
tering time t as well as short mean free path I = y/ Dt, 
then the results are not sensitive to I and r, which are 
fixed throughout the paper. Note that similar crossover 
also happened in graphene but is completely opposite 
[26l [34] because graphene is gapless and has ignorable 
spin-orbit interaction [35] 136] . 

The temperature dependence of the quantum conduc- 
tivity a F also provides the signature of the crossover. Be- 
cause the weak (anti-)localization comes from the quan- 
tum interference, their temperature dependence mainly 
originates from the phase coherence length Ia. In the 
quantum diffusion regime, becomes much longer than 
the mean free path i and increases with decreasing 
temperature as decoherence mechanisms from electron- 
phonon and electron-electron interactions are suppressed. 
The similarity of band dispersion to the conventional 
two-dimensional systems [37] allows us to assume that 
U = C ph {T/T*)-V 2 + CeeiT/T*)- 1 / 2 , where the coeffi- 
cients C p h from electron-phonon and C ee from electron- 
electron interactions are determined by imposing that 
£ = 300 nm at T/T* = 0.01 and 1$ = I at T/T* = I. T* 
is a characteristic temperature above which the system 
transits from the quantum to classical diffusion regimes, 
so a F vanishes at T = T* . The order of magnitude for 
T* ranges between 10 and 100 K [37H39] . Figure [3] shows 
<J F (T) in Eq. Q for different tq/t\. The crossover from 
WL to WAL can be read from the sign of a F turning from 
negative to positive. Note that the total conductivity is 
a(T) = a D (T)+a F (T, B), where the Drude conductivity 
a D may give extra temperature power law that depends 
on scattering mechanisms 38, 39 . But a F from the weak 
(anti-)localization is the contribution most sensitive to 
magnetic field (because of the logarithmic dependence). 
<j F thus can be extracted from the total conductivity by 
a F (T) = a(T,B c ) - a(T,0), with B c a finite magnetic 
field large enough to quench the weak (anti-)localization, 
say, of 1 Tesla. 

We can also have some qualitative arguments for the 



conduction bands. All four conduction bands can take 
part in transport as they always intersect the Fermi sur- 
face together. Therefore, different from the transport 
in hole-doped samples, both intravalley spin-orbit and 
intervalley spin-conserved scattering are possible in the 
conduction bands. The rule is that spin-flip scattering as- 
sisted by spin-orbit interaction leads to WAL while spin- 
conserved scattering to WL because the pseudo-spin de- 
grees of freedom is frozen. All types of scattering and the 
resulting weak (anti-)localization are shown in Fig. [4] 
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